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ABSTRACT 

l/~) ' Context. Classical T Tauri stars (CTTS) differ in their X-ray signatures from older pre-main sequence stars, e.g. weak- 

lined TTS (WTTS). CTTS show a soft excess and deviations from the low-density coronal limit in the He-like triplets. 
Aims. We test whether these features correlate with accretion or the presence of a disk by observing IM Lup, a disk- 

i i bearing object apparently in transition between CTTS and WTTS without obvious accretion. 

("V | Methods. We analyse a Chandra grating spectrum and additional XMM-Newton data of IM Lup and accompanying 

X_/\ , optical spectra, some of them taken simultaneously to the X-ray observations. We fit the X-ray emission lines and 

decompose the Ha emission line in different components. 

i-C ■ Results. In X-rays IM Lup has a bright and hot active corona, where elements of low first-ionisation potential are 

ft' depleted. The He- like Neix triplet is in the low-density state, but due to the small number of counts a high-density 

C*. [ scenario cannot be excluded on the 90% confidence level. In all X-ray properties IM Lup resembles a main-sequence 

star, but it is also compatible with CTTS signatures on the 90% confidence level, thus we cannot decide if the soft 
excess and deviations from the low-density coronal limit in the He-like triplets in CTTS require accretion or only the 
presence of a disk. IM Lup is chromospherically active, which explains most of the emission in Ha. Despite its low 
equivalent width, the complexity of the Ha line profile is reminiscent of CTTS. We present an estimate for the mass 
accretion rate of 10 _11 Mq yr _1 . 
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on: 

1. Introduction ()Fang et alJ feOO!)). Usually the combination of strong Ha 

' . . . . . emi ssion and line asymmetry is a reliable accretion indica- 

• . Our view on the formation of stars and planetary systems tQr (M uzcrollTJaS][l998, 200J). 

c and their emergence from molecular clouds has made sigmf- WTTg with thdr 1qw Ha EW &re all expecte d 

icant progress over the last decades. Stars form when giant nQt tQ show accretion . It is sugge stive to interpret the 

_H , molecular clouds fragment and contract to form proto-stars. WTTg ag more evolyed CTTgj where accretion has stopped 

.. . Mass accretion onto those stellar cores proceeds via an ac- already This doeg nQt necessarily impl the absence of a 

. £ ; cretion disk while the surrounding envelope eventually dis- d[sK ag transitional disks may exist without ongoing ac- 

K> . perses. The low-mass pre-mam sequence stars in this s stage cretion for SQme time (Padgctt ^^M) ■ For many TTS 

H ■ c ° me m ^° ^^M h e classical T Tauri s tars (CTTS) and the Ra EW ig knQwn tQ be variabl6; possibl becauge the 

_ the weak-lined T Tauri stars (WTTS). Traditionally, they accretion switches on and off or _ at least _ the accre tion 

were distinguished only by their Ha emission line equiv- , i 

? 6n * W ^c th ; Wlth t ?° Se St f S °I E ^^ l° u A d6fined t0 , ' Both types of TTS have been known for a long t ime to 
be CTTS It turned out that the EW of Ha is a good be c ioug x _ r emitterg d FeigeLson fc MontmerleHlQQSh . 

™c ^ acCTetl Z^°Y T b °^ h l Z^l u ! As a class CTTS stand out fro m other X-ray sources by 

CTTS still possess. The Ha EW of WTTS is smaller and their strong soft X -ray excess (|Robrade fc SchmittJ 1200% 

their line profiles are symmetric; in contrast CTTS exhibit ^del^Tcllcschi 2007), but CTTS also exhibit normal 

broader emission lines, which are sometimes asymmetric CQronal activity &nd gteUar flar ^ 

Send offprint reauests to: H. M. Giinther, T^SFP ^ S^ ^^ llonn J iF^™^^^ 

e-mail: moritz . guentherOhs . uni-hamburg . de and XMM-Newton (|Kastner et al.||2002|; IStelzer fc Schmittj 

* Based on observations obtained with XMM-Newton, an ES A 12QM), was the first star where unusual line ratios in the 

science mission, and Chandra, a NASA science mission, both He-like triplets was found indicating high densities in the 

with instruments and contributions directly funded by ESA formation region. This phen omenon repeats in most CTTS 

Member States and NASA . observed so far (BP Tau: ISchmitt et aD ()2005D : V4046 



CZ2 
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See: IGiinther et all (|2006h: RU Lup: iRobrade fc Schmittl 
(l2007h: MP Mus: lArgiroffi et all fl2007Tl : Hen 3-600: 
Huc nemoerder et all (|2007l )). There are exceptions to this 
rule - in the more massive, eponymous T Tau itself, 
although known for its high accretion rat e, the O vh 
triplet is consistent with the coronal limit (|Giidel et all 
l2007bh . The same has been found for HAcBe stars, which 
are in a similar evol utionary state as the C TTS, yet of 
higher mass (AB Aur: iTelleschi et all (|2007cD : HD 163296: 
IGiinther fc Schmittl (|2009ft l The soft excess and the He- 
like triplet line ratios in CTTS can be natu rally linked to 
mass accretion from proto-planetary disks (|Lamzinlll998t 
IGiinther et all l2007h . 

Systems containi ng typical WT TS such as TWA 4 
(|Kastner et all l200l and TWA 5 (|Argiroffi et all [2005h 
show no sign of active accretion and no high densities in 
their X-ray spectra. Well studied examples of young main- 
sequence (MS}_stars_such as AU Mic, Speedy Mic and 
AB Dor ( Nes s et alll2004h also show triplets indicating low 
densities sim ilar to older M S obje cts. Still, 5 of 83 WTTS 
observed by iPadgett et all (|2006l ) do show an IR-excess 
caused by dust within a few AU from the star, detected 
by Spitzer, thus there is a matter reservoir for accretion. 

Observationally, we could not distinguish in the past 
if the X-ray signatures observed in CTTS in contrast to 
WTTS and young MS stars are caused by active accre- 
tion or the presence of a disk. To solve this question we 
observed a system with Chandra, that was previously clas- 
sified as disk-bearing, non-accreting WTTS. Simultaneous 
to the X-ray observations, the accretion state is observed 
using optical spectroscopy. 

The properties of the our target IM Lup are sum- 
marised in Sect. [5] We then present the X-ray observations 
of IM Lup and the accompanying optical data in Sect. [3] 
We show our new results and compare IM Lup to CTTS, 
WTTS and MS stars in Sect. |U In Sect. Owe discuss the 
implications and end with a short summary in Sect. [6] 

2. Properties of IM Lup 

In this article we present new observations of IM Lup 
(=Sz 82), the only known X-ray bright transition object 
with a low Ha EW (4 A) (|Wichmann et al.l[l999h . far less 
than the values for typical CTTS such as BP Tau and 
TW Hya. Still, the IR-excess in IM Lup indicates the pres- 
ence of a disk (jPadeett et alJl20 06jl making it a favourite 
object for disk evo l ution studies in the I R (starting with 
IPadgett et ail 120061: ISchegerer et alll2006D . 

There is good evidence that IM Lup is a s ingle object 
both from ground-based speckle observ ations (iGhez et all 
19971) and from space-b ased HST data (jPinte et alj |2008). 



Table 1. Observing log of X-ray observations. 



Observatory 


ObsID 


Obs. date 


Exp. time 


XMM-Newton 


0303900301 


2005-08-08 


30 ks 


XMM-Newton 


0303900401 


2005-08-17 


27 ks 


XMM-Newton 


0303900501 


2005-09-06 


27 ks 


Chandra 


9938 


2009-05-04 


150 ks 



Wic hmann et all fl998) derived a HIPPARCOS distance of 



190 ± 27 pc, larger than earlier estimates, but still consis- 
tent to them within the errors. We use this value through- 
out the paper. The spectral type of the star is M0, how- 
ever, there is some debate on its mass. The values vary 
with the pre-main se quence evolutionary tracks adopted: 
Based on the data of Hughes et all (|l994fl the models of 
iD'Antona fc Mazzitellil (ll994ri yield M = 0.29 M^ and the 
models of lSiess et all (|2000D give M = 0.52M^.lPinte et all 
(|2008l) use the tracks from fearaffe et all (| 19981 ) and find 
M = IM© and R = 3R@. The age is around 1 Myr in 
all models. Observationally there is support for a mass of 
1.2 ± 0AM(7) from millimetre observations, which 



M 



'0 



resolve the emission lines from a disk in Ke plerian rota- 
tion; the inclination of this disk is 54° ± 3° (|Panic et all 
2009 ). The redden ing has been derived between Ay = 0.5 
(jPinte et alll2008l) and A v = l (|Hughes et al.lll994h . 

The inner disk of IM Lup is resolved down to < 0.4 AU, 
dust grains have partially settled in the disk-midplane and 
grown to a few millimetres; the disk surface shows strong 
silicate emission bands and there are sig natures of a low 
degree of crystallisation (|Pinte et alll2008l ) . Beyond 400 AU 
the surface density of CO gas and dust grains with sizes 
1 Mm- lm m drops dramatica lly. It remains detectable out to 
900 AU (panic et alJl2009h . 



3. Observations and data reduction 

3.1. X-ray observations 

IM Lup was observed off-axis three times with XMM- 
Newton giving a total exposure time of 84 ks and with 
Chandra for 150 ks in a dedicated observation using the 
high-energy transmission grating (HETGS). The observa- 
tion log is summarised in table [TJ All data was reduced 
using the standard XMM-Newton Science Analysis System 
(SAS) software, version 9.0, (jGabriel et all .2004) or the 
Chandra Interactive Analysis of Observatio ns (CIAO) soft- 
ware, version 4.1.2 (|Fruscione et all 120061) . in the case of 
the Chandra observation. The XMM-Newton observations 
are contaminated by proton flares and we applied the stan- 
dard selection criteria. We extracted CCD spectra for the 
PN and the MOS2 detector. During the observations one 
of the chips in the MOS1 camera was not operating and, 
unfortunately, IM Lup is located in this gap for two of the 
three observations. We primarily used the MOS2 data for 
lightcurves because it is less affected by the proton noise 
and the PN data for spectra because of their higher signal. 
We extracted lightcurves and CCD spectra from all three 
observations. We did not reprocess the Chandra evt2 file, 
thus the predefined selection criteria apply. Positive and 
negative orders for the HEG and MEG are merged. Spectral 
fittin g was carried out using XSPEC V12.5.0 (jArnaudl 
119961 ). and individual li ne fluxes were measured u sing the 
CORA line fitting tool (|Ness fc WichmaniJl2002h . CORA 
employs a maximum likelihood method, taking the Poisson 
statistics into account. Because the line widths are domi- 
nated by instrumental broadening, we keep them fixed at 
AA = 0.02 A for the MEG and AA = 0.01 A for the HEG, 
the line profile is given by a modified Lorentz profile with 
j3 = 2.5. We fix the wavelengths at the theoretical values 
taken from the C HIANTI 5.1 database (jDere et all 119971 
lLandi et all 120 06!), but by fitting the strongest lines with 
a free wavelength we verified that our wavelength scale is 
accurate. The error introduced by fixing the wavelength 
is negligible compared to the uncertainties from counting 
statistics in all cases. 
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Table 2. Log of optical observations. 



No 


Observatory 


ID 


Obs. time 


Exp. time 


HI 


ESO/HARPS 


6201 


2008-05-09 02:19 


20 min 


H2 


ESO/HARPS 


6381 


2008-05-09 08:55 


20 min 


H3 


ESO/HAPRS 


6711 


2008-05-10 02:18 


20 min 


H4 


ESO/HARPS 


6971 


2008-05-10 08:10 


20 min 


Al 


ANU 


35 


2009-05-05 12:38 


30 min 


A2 


ANU 


36 


2009-05-05 13:09 


30 min 


A3 


ANU 


36 


2009-05-06 11:50 


30 min 


A4 


ANU 


37 


2009-05-06 12:20 


30 min 



3.2. Optical observations 

Observations were carried out using the 79g/mm echelle 
spectrograph (R« 24, 000) on the Nasmyth-B focus of the 
Australian National University's 2.3m telescope (ANU) on 
the nights of May 5th and 6th, 2009. Each night two 1800 s 
exposures were taken one after another. The observing log 
for all optical observations is shown in table [2 A GG490 fil- 
ter is used in the beam to avoid overlapping of the first two 
orders at the reddest wavelength. A cross disperser grat- 
ing of 316/750nm was used to allow collection of data of 
orders 55-34 or 4120-6630 A. Bias frames, ThAr Arc for 
wavelength calibration and quartz lamp exposures for flat- 
fielding were also included for the process of data reduc- 
tion. The code used to reduced the 2.3m data is a slightly 
modified version o f ESpRIT (Echelle Spectra Reduction: an 
Interactive Tool) (jDonati et al.lll997| ). 

For comparison purposes we also downloaded observa- 
tions from the ESO 3.6m telescope at La Silla obtained 
with HARPS for IM Lup from the ESO archive (program 
ID 081.C-0779(A)). These datasets are fully reduced. The 
uncertainty of the flux in these observations was taken as 
the Poisson error of the count statistic. This is consistent 
with an estimate from the spectra, where we chose an ap- 
parently line-free region close to the Ha line. We calculate 
the mean flux and the standard deviation. The relative er- 
ror in the spectra is 3% for each bin in the oversampled 
spectrum, that is a S/N about 100 rebinned to instrumen- 
tal resolution, except for the second exposure, where the 
relative error is 5%. 

To compare IM Lup to MS stars, we obtained spec- 
tra of HP 156274 and HD 156026 from the library of 
iBagnuloetall (120031). 



4. Results 

4.1. The X-ray lightcurve 

In Fig. [1] the X-ray lightcurve of IM Lup is shown. The 
XMM-Newton observations are separated by several days 
and IM Lup shows a steady increase in luminosity during 
this time, holding the ratios between hard (1.5-10.0 keV) 
and soft (0.3-1.5 keV) emission nearly constant. Due to dif- 
ferent detector efficiencies the hardness ratios of Chandra 
and XMM-Newton cannot be compared directly. During 
the much longer Chandra observation, there is small scale 
variability in both bands, but without clear flaring activity. 
A KS-test shows, that a constant lightcurve can be rejected 
on the 72 % level for the soft band and on the 98 % level 
for the hard band. 
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Fig. 1. Lightcurve in soft (black) and hard (red/grey) band. 
Due to different detector efficiencies XMM-Newton and 
Chandra cannot be compared directly. The shaded areas 
show the times of the optical observations. 



0.2 




Energy (keV) 



Fig. 2. CCD spectra of IM Lup from XMM-Newton 
EPIC/PN (ObsID 0303900401, upper curve) and Chandra 
ACIS (ObsID 9938, lower curve) with the best-fit model 
overlaid. See table |3] for model parameters. 



4.2. The CCD X-ray spectra 

In Fig.[5]we show the MOS2 CCD spectrum from the second 
XMM-Newton exposure and the zeroth order ACIS spec- 
trum, both binned to a minimum of 25 counts per bin. At 
first sight we can already identify the Nex Lya line, caus- 
ing a peak around 1 keV and the Sixni and Sixiv emis- 
sion around 2 keV. We describe the plasma in more detail 
with three optically thin emission components according 
to the vapec model combined with a single photoabsorp- 
tion model in XSPEC and overplot the best-fit model. A 
simpler model with only two components and coupled ab- 
sorbing column density njj, temperature, and abundances 
for all four observations already provides a statistically ac- 
ceptable fit (red. x 2 = 1.12). However, the residuals show 
systematic trends for low energies indicating an inaccu- 
rate combination of nn and T\. Thus we added a third 
vapec component, keeping the temperatures of the three 
components coupled over all datasets. This ensures that 
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the relative contribution of soft, medium and hard compo- 
nent can be compared. There is no significant difference in 
the fitted njj between the separate datasets, thus we also 
keep it coupled. The parameters of our best-fit model are 
given in table [3J The errors are calculated as 90 % con- 
fidence ranges and t he fitted abundances are relative to 
the solar values from iGrevesse fe Sauvall (J1998D : all other 
abundances are fixed at their respective solar value. Some 
of the values in the table are strongly coupled. In order 
to retain the overall luminosity of the star, when e.g. the 
soft component is lowered the emission measure (EM) of 
the medium component has to be increased. In all observa- 
tions the emission is dominated by a hot plasma component 
around kT = 2.1 keV. The separation between the consecu- 
tive XMM-Newton exposures is 9 and 19 days. In this time 
the total plasma emission measure increases by 60%; at the 
same time a larger portion of the total emission measure 
is concentrated at higher temperatures. The lightcurves in 
Fig. [T] show, that this increase is steady during the obser- 
vations. In the Chandra observation the distribution of the 
emission measure in the different components is comparable 
to the second XMM-Newton dataset, only the hot compo- 
nent is reduced by about a third; thus, the total emission 
measure is lower than in any of the XMM-Newton obser- 
vations. We measure the absorbing column density towards 



Table 3. Best-fit model for the CCD spectra (90% conf. 
ranges) 



ObsID 



XMM-Newton Chandra 

...301 ...401 ...501 9938 



Abundances 



O 

Ne 

Mg 

Si 

S 

Fc 



nk [10 2i cm" 2 ] 
kTi [keV] 
kT 2 [keV] 
kT 3 [keV] 
EM X [10 52 cm" 
EM 2 [10 52 cm" 
EM 3 hn52 



-+U.3 
-0.2 

2.o±g:« 

0.8+°;* 
0.2±g;| 

1+0.2 

^-O.l 

+0.2 

0.1 



0.5 4 





0.4: 



1.3 ±0.2 
3+ 01 

u -°-0.3 



0.6 
2.1 



[10 5 



cm 



6i 9 6 

5l 8 5 
34« 



3 
10 

38 



13 

3 

8 

10 

+6 

5 



+0.1 

0.1 
+0.2 
-0.1 

Q+35 

12±* 



60 



+7 

7 



log L x [erg s 

T2 



unabsorbed luminosities (0.3-10 keV) 

TT 



o+8 
°-3 

ioiI 

+4 

3 



26 



30.7 



30.8 



30.9 



30.7 



X A = 266 for 254 dof (red. \ l = 1-05) 



the source as A^h = 1.3 ± 0.2 x 10 21 cm -3 . We use this 
value to calculate the total unabsorbed X-ray luminosities 
Lx given in table |3j Assuming a standard gas-to-dust ra- 
tio the optical reddening and the X-ray absorption should 
be related through the formula A^h = Ay ■ 2 x 10 21 cm~ 2 
(jVuong et al.l 120031 ) . leading to Ay = 0.7 mag in a gree- 
ment with t he va lues found by iHughes et al.l (|1994h and 
iPinte et al.l (|2008t ). 



4.3. The high-resolution X-ray spectra 

4.3.1. Line fluxes 

The MEG grating spectrum is shown in Fig. [3] We mea- 
sure line fluxes from the HEG and the MEG. In table [4] 
we give fluxes for all lines which were detected above the 
2<j level. In the case of the He-like triplets we quote all 
lines, if one of them is strong enough to be included in 
the table. From the fitted line fluxes we calculate the to- 
tal line luminosity, using the absorption cross section from 
Balucinska-Church & Mc Cammonl (jl992f ) . The relative er- 
rors on the line intensity are larger than those on the line 
fluxes, because here the uncertainty in the distance and the 
fitted Ah contribute to the total error budget. 

Table 4. Measured line fluxes (la errors) 



Ion 



"MAT 



counts photon flux" luminosity 



HEG 



Si xiv 
Mgxn 

Nex 
Nex 
Nex 



6.18 

8.42 

9.71 

10.24 

12.13 



7.6 ± 3.3 

9.0 ±3.4 

6.1 ±2.9 
5.6 ±2.5 

30.5 ±5.7 



1.2 
1.5 
1.8 

2.6 
20.3 



±0.5 
±0.6 
±0.8 
±1.2 
±3.8 



17 ±9 

17 ±8 

19 ±10 

27 ±14 

197 ± 68 



MEG 



Si xiv 

Si xiii 

Si xiii 

Si xiii 

Mgxn 

Nex 

Nex 

Fe xviii 

Fe xvi 1 1 

Nex 

Fe xvn 

Neix 

Neix 

Neix 

Ovin 



6.18 

6.65 

6.69 

6.74 

8.42 

9.71 

10.24 

11.25 

11.53 

12.13 

12.26 

13.45 

13.56 

13.70 

18.97 



12.2 ±4.7 
18.0 ±5.2 

7.5 ±4.0 

15.3 ±5.1 
14.9 ±5.0 
14.0 ±4.7 
22.6 ±5.7 

6.7 ±3.0 
8.9 ±3.4 

90.8 ±9.7 
10.2 ±3.8 

20.9 ± 5.0 
3.3 ±2.6 

10.9 ±3.9 
7.9 ±3.3 



0.9 

1.2 
0.5 
0.9 
0.9 
1.4 
2.6 
1.1 
1.5 
19.7 
2.3 
7.1 
1.2 
4.0 
9.9 



±0.3 

±0.4 
±0.3 
±0.3 
±0.3 
±0.5 
±0.6 
±0.5 
±0.6 
±2.1 
±0.9 
±1.7 
±0.9 
±1.4 
±4.1 



13 ±6 

17±7 

7±4 

12 ±5 

10 ±4 

14 ±6 
26 ±10 

11 ±6 

15 ±7 
191 ± 59 

22 ±10 
70 ±26 
12 ±10 
40 ±18 
103 ± 54 



a in units of 10~ 6 s" 1 cm" 2 

b unabsorbed, in units of 10 27 erg s _1 

We detect lines of O, Ne, Mg, Si and Fe: The first four 
elements turn up in the Lyman series and in the He-like 
triplets; iron can be seen in the ionisation stages Fexvn 
and Fe xviii. Due to the low signal only five lines can be 
detected in the HEG and all fluxes agree with the MEG 
measurement within the statistical uncertainty. In the fol- 
lowing we use only the MEG data. 

A special analysis is required for the Ovu He-like 
triplet, because its signal is very weak due to absorption 
and Chandra's low effective area in this region. Its flux 
is important to compare IM Lup to other CTTS and MS 
stars. There are six photons recorded between 21.5 A and 
22.2 A (at 21.60, 21.72, 21.80, 22.06, 22.07 and 22.10 A), 
which is compatible with the background of 13 cts A -1 
determined from the line-free region from 20 A to 23 A 
around the O vii He-like triplet. However, for the line pro- 
file of Chandra 88% of all photons are expected within one 
FWHM=0.02 A. Three of the six photons fall within this 
range of the theoretical wavelengths (21.6020, 21.8071 and 
22.1012 A according to the CHIANTI database), but only 
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Fig. 3. MEG spectrum of IM Lup, rebinned by 2. The lines detected with more than 2a significance are marked (table0]). 
The inset shows an enlargement of the He- like Neix triplet. 



0.25 photons are expected from the background in the cen- 
tral FWHM of each line. According to Poisson statistics 3 
photons on a background of 0.75 photons are a detection of 
O VII emission at the 95% significance level, but the uncer- 
tainty on the derived flux is obviously very large. The 90% 
upper limit on the flux, i.e. the flux, where the chance is 
10% to observe three or less photons in the O vn triplet is 
2 x 10~ 5 cts s _1 cm~ 2 , which corresponds to a dereddened 
luminosity of 3 x 10 29 erg s~ 1 . 



4.3.2. He-like triplets 

He- like ions show a triplet of lines, where the line ratios are 
temperature and density sensitive. These triplets consist 



of a resonance 



an intcrcombination (i) and a forbid- 



den line (/) (jGabriel fc Jordan! 1 19691: iPorquet et al.ll2001h . 
Commonly the R- and G-ratios (R = f/i and G = {f+i)/r) 
are used; for high electron densities n c or strong UV photon 
fields the i?-ratio falls below its low-density limit, because 
electrons are collisionally or radiatively excited from the up- 
per level of the / to the i line, but the UV field of late-type 
stars like IM Lup is far too weak to influence the i?-ratio. 
We interpret the ratios with the CHIANTI 5.1 database. 

According to table 2] the observed i?-ratio for Sixm is 
1.8, with a 90 % lower limit of 1.0 according to a Monte- 
Carlo simulation of Poisson-distributed counts using the 
obs erved count numbers as e xpectation values (for details 
see I Gunther fc Schmitd I2009L appendix A). The predicted 
low-density limit for this ratio is 2.0 with smaller num- 
bers for densities logn c < 12.5. Not surprisingly the Sixm 
triplet is compatible with the low density limit (90 % range: 
logn c < 13.8). The G-ratio is temperature sensitive. For 
temperatures, where more than 10% of the Si is present 
in the form of Sixm, its value is predicted between 1.2 
(logT = 6.1) and 0.6 (logT = 7.2). This all falls in the la 
error range of the observations, thus the G-ratio does not 
constrain the models. 

In the Neix triplet the formal //z-ratio is 3, the Monte- 
Carlo simulation gives lower boundaries on the ratio of 2 
(la) and 0.8 (90%). According to CHIANTI the f/i ratio in 
the low density regime is 2.8, thus the respective limits on 



the electron density in the emitting region are logn c < 11.4 
(la) and logn c < 12.3 (90%). The intercombination line in 
the Ne IX triplet is known to be contaminated with Fe lines 
in many cases, however, iron is under abundant in IM Lup 
and this would further increase the f/i ratio. The Neix r 
line also contains a blend from Fe xix. We cannot correct for 
this blend because no single Fexix is detected and, given 
the low abundance of Fe, we estimate that this blend is far 
smaller than the observational uncertainty of the line flux. 
In summary, the R-ratio of the Si xm triplet limits the 
electron density of the emitting plasma to n c < 13.8, the 
R-ratio of the Neix triplet to logn c < 12.3 (90% confidence 
in both cases). 



4.4. The optical spectra 

The optical spectra are dominated by photospheric ab- 
sorption lines. From the high-resolution HARPS data we 
measure the rotational broadening to vsini = 15 ± 
2 km s -1 by comparison t o synthetic PHOENIX spectra 
(jHauschildt fc Bar on 1999). Using the st e llar r adius and 
the inclination determined bv iPinte et a l. (2008) this sug- 
gests a stellar rotation period of 8 days if the disk incli- 
nation and the stellar inclination match. We compared the 
spectra to two templates taken from MS stars of similar 
sp ectral type (HD 156274 and HD 156026) in the library 
of iBagnulo et al.l (2003). The spectra are broadened and 
degraded to the spectral resolution of our IM Lup data. 
The depth of most photospheric lines matches well, thus 
no significant veiling is present in IM Lup. 

In the following we discuss where IM Lup differs from 
the templates. We find Lil in absorption at 6708 A, con- 
firming its youth. The equivalent width of this line is about 
0.3 A. Typical activity tracers like Ca K and the Balmer 
series are in emission. Figure S] shows the Ca H and K lines 
and Fig. [5] the Ha line profile in the four different expo- 
sures taken simultaneously to the Chandra observation. 
The Ha equivalent width increases during our observations 
from 2.8 A on the first day (HI and H2) to 4.5 (H3) and 
5.2 A (H4) in the exposures on the second day, respectively, 
where we measured the emission core of the line only (in 
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Fig. 4. Ca H and K lines in all four exposures. For both 
lines only sharp emission peaks are visible. The small line 
at 3970 A is He. 




2009-05-05 12:38 

i i I i i 

-500 



500 



ine shift [km/sl 



Fig. 5. Ha line profile in all four exposures. The line is in 
emission and asymmetric. A broad absorption component is 
also visible for the second day. Observations are continuum 
normalised and shifted in intensity for clarity. The red/grey 
data points are excluded from the fit. Solid lines show a fit. 
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Fig. 6. Ha line profile in four HARPS exposures one year 
before our observations. Observations are continuum nor- 
malised and shifted in intensity for clarity. The red/grey 
data points are excluded from the fit. The lines show the 
best fit of five Gaussians according to table [SJ The single 
components are indicated for the first observation. Note the 
smaller velocity scale compared to Fig. \5\ 



accordance with the X-ray data we use the convention to 
denote emission lines with a positive EW and absorption 
lines with negative values). 

In all cases the Ha EW stays below 10 A, which is the 
canonical dividing line between CTTS and WTTS. The line 
profile is asymmetric, with additional emission on the blue 
side. On the first day a small bump exists at +100 km s -1 , 
which might be due to an unresolved secondary peak at 
this velocity. This seems to be t he typical line pr o file fo r 
IM Lup; it was also observed bv IWichmann et al.l ([1999). 
In a HARPS observation (Fig. [6]) secondary peaks on both 
sides of the Ha line appear. On the second day of our 
ANU observations (A3 and A4) broad absorption wings are 
present, which are only marginally visible on the first day 
(Al and A2), while at the same time the central emission 
is strong er peaked. Acc ording to the classification scheme 
oflReipurth et alJ (Il996h the Ha line profile in 1998 is III-R 
(|Wichmann et al.lll999h . in 2008 III-B/R, in 2009 III-R on 
the first day and IV-B/R on the second day. 
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In parallel to Ha the other detected Balmcr lines and 
also He I 5867 A increase in strength on the second day in 
contrast to the X-ray emission. Surprisingly, Fig. [T] shows 
that the X-ray count rate during the observations on the 
second day is only half of that on the first day; the hardness 
ratio is constant. 

For comparison we analyse a HARPS dataset of four 
exposures of IM Lup, which was taken one year before 
the Chandra observations. The Ha emission line, shown in 
Fig. |51 looks similar; its total equivalent width (including 
the absorption component) decreases during the observa- 
tions from 4.3 A to 2.7 A. The spectral resolution in these 
observations is better, thus we can fit the separate compo- 
nents of the line reliably. We exclude two narrow absorp- 
tion lines from the profile fitting (Til and Fei) and also 
ignore the center of the Ha line, which is likely affected 
by the optical depth. We find that in total five Gaussian 
components are necessary to describe the line profile: The 
most prominent component is the central emission peak, 
then there are two narrow emission components shifted by 
about 100 km s _1 to the blue and red side, respectively, 
and two broad components, one of them in absorption. 
The components are listed in table [5] and shown in Fig. [5] 
Absorption below the continuum is only seen for velocities 
above 140 km s _1 . Still, the wings on the red side are best 
fit by a wide absorption component with a centroid close to 
the stellar rest velocity, if we prescribe a Gaussian shape. 
If the absorption on the red side in non-Gaussian, than the 
wide emission component is no longer needed and there is 
no absorption on the blue side of the line, so that the three 
narrow emission components are sufficient to describe the 
line shape. The central emission matches the stellar rest ve- 
locity at all times, its amplitude changes with time. In the 
fit this is seen as decreasing flux in the wide component, 
while the central component is nearly constant. The narrow 
red-shifted component moves towards the center during the 
observations, while the narrow, blue-shifted emission keeps 
its position. 

In the ANU observations a qualitatively similar struc- 
ture is observed, but the resolution and signal-to-noise ra- 
tio is lower, so only two Gaussians are required. Absorption 
seen on both sides of the line on the second day of the ANU 
observations (A3 and A4 in Fig. [5]) is consistent with a wide 
absorption component centred on the stellar rest velocity 
and sigificantly broader than in the HARPS dataset. The 
FWHM of the main emission feature of 160 km s _1 and 
the velocity shift of the bump on the red side of the line 
profile in the first two observations agrees with the nar- 
row red component in the HARPS observations (table [5]). 
The absorption in the ANU observations is stronger and 
broader than in the HARPS data, in fact, it is stronger 
than the emission line on the second day. Gaussian profiles 
are not ideal to fit the line in Fig. \E\ the fit is too shallow 
and too broad on the red side around 150 km s _1 , but they 
illustrate the change in relative strength between the ob- 
servations. The fits underpredict the flux in the line center, 
because the line is more strongly peaked than a Gaussian. 



4.5. Comparison X-ray properties to CTTS, WTTS and MS 
stars 

To compare the X-ray luminosity with other TTS of similar 
spectral type we used the data published from the COUP 



Table 5. Properties of line profile components 



Comp. 


No 


EW° 


centre 


FWHM 






[A] 


[km s^ 1 ] 


[km s- 1 ] 


HARPS - fit with 5 Gaussians 


central 


HI 


3.9 ±0.1 


-2.6 ±0.2 


76 ±1 




H2 


3.6 ±0.2 


-1.7 ±0.4 


77 ±2 




H3 


3.9 ±0.1 


-5.7 ±0.2 


85 ±1 




H4 


3.9 ±0.1 


-4.2 ±0.3 


82 ±1 


blue wide 


HI 


3.4 ±1.1 


-62 ±9 


248 ± 13 




H2 


4.4 ±3.8 


-53 ± 22 


249 ± 30 




H3 


0.8 ±0.1 


-133 ±3 


115 ±4 




H4 


1.5 ±0.4 


-128 ±7 


166 ± 10 


blue narrow 


HI 


0.34 ± 0.02 


-92.2 ± 0.6 


51 ±2 




H2 


0.29 ± 0.04 


-90.5 ±1.2 


43 ±4 




H3 


0.55 ± 0.02 


-95.6 ± 0.4 


49 ±1 




H4 


0.62 ± 0.04 


-95.5 ±0.5 


54 ±2 


wide abs 


HI 


-3.9 ±2.3 


-12±4 


412 ±24 




H2 


-4.6 ±5.3 


-5 ±16 


354 ± 43 




H3 


-2.2 ±0.3 


36 ±3 


312 ±6 




H4 


-3.8 ±0.7 


7±4 


351 ±7 


red narrow 


HI 


0.61 ±0.04 


100.6 ± 0.5 


60 ±2 




H2 


0.51 ±0.09 


96.3 ±1.4 


60 ±5 




H3 


0.48 ±0.01 


93.8 ±0.4 


55 ±1 




H4 


0.47 ± 0.03 


86.2 ±0.9 


62 ±2 




ANU 


- fit with 1 or 


2 Gaussians 




central 


Al 


2.8 ±0.2 


-25.8 ±1.0 


159. ± 3. 




A2 


2.8 ±0.2 


-27.5 ±1.0 


157. ±2. 




A3 


4.5 ±0.2 


-22.6 ± 0.6 


161. ±2. 




A4 


5.2 ±0.2 


-21.7 ±0.5 


150. ±1. 


absorption 


A3 


-4.7 ±2.8 


15. ±7. 


750. ± 50. 




A4 


-7.8 ±5.1 


-36. ±6. 


1050. ± 50. 



absorption is shown with negative values 



dPreibisch et all [2001 and the XEST (iGiidel et al.ll2007af) 
projects. Both are large X-ray surveys covering the star 
forming regions of Orion and the Taurus molecular cloud, 
respectively. IM Lup is a bright object compared to both 
datasets. 

Compared to the Taurus molecular cloud 
(JTelleschi et al.l l2007ah the mean temperature and X- 
ray luminosity Lx of IM Lup would make it the most 
luminous CTTS. If it is a WTTS it still belongs to the 
brightest 20% of all WTTS in that region. 

The bolqmetri c luminosity of IM Lup is L& ; = 1.9Lq 
(|Pinte et al.ll2008b so log(L x /L bol ) = -3.1. This number is 
typical both for MS stars, which satu rate at log(L j y /Lboi) = 
-3, and also for CT TS or WTTS dFlaccomio et al.1 12003: 
lPreibischet~aTll2005ll . 

The abundances of IM Lup shown in table [3] exhibit an 
IFIP (inverse first ionisation potential) pattern, where ele- 
ments of large first ionisation potential are enriched relative 
to the solar abundance and those of low FIP are depleted. 
This is observed in most active stars, independent of their 
evolutionary state. The values found in IM Lup agree wit h 
typical abundances as given in the review by Giidcl (2004). 

iRobrade fc Schmittl (l2007ll and IGiidel fc Telleschil 
(2007) present an analysis of the O Vlll/O VII ratio as a 
measure of the excess of soft emission in CTTS with re- 
spect to stars on the MS. These studies show the observed 
excess to be confined to a narrow temperature range at 
the formation of the He-like O vn triplet at 1-2 MK. The 
upper limit we determine for the O vn emission in IM Lup 
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is compatible with both, MS stars and the soft excess in 
CTTS. 

The f/i ratio of the Neix triplet in IM Lup is sta- 
tistically compatible on the la level with typical values 
of coronal sources (jNess et al.l 120041 ) in the low-density 
limit; however, a reduced f/i value is allowed for IM Lup 
on the 90% confidence level due to the low count statis- 
tic. The l ow-density limit is also obser ved in WTTS 
(TWA 4: (jKastner et all J200l; TWA 5: (lArgiroffi et al.l 
[2005h and HDE 245059 dBaldovin-Saavedra et all 120091) ). 
Most CTTS on the other hand show a reduced f/i ratio 
which is at tributed to the i r accretion spots on the surface . 
TW Hva dKastner et all 12003 iStelzer fc Scimlittl l200l 
BPTau (ISchmitt et all 120051). V4046 See (iGunther et al.1 
[200l. RU Lup ([ Robradefe Schmittl I2007T). MP Mus 



(lArgiroffi et al.ll2007t) and Hen 3-600 ()Huenemoerder et al.l 
2007) all show significant deviations from the low-density 
limit. 

Going beyond the He-like triplets we compared our 
MEG spectrum (Fig. H) o f IM Lup with a Chandra/MEG 
spectrum of YY Gem (JLiefke et al.l [200811 . a MS binary of 
spectral type M3-M4, similar to IM Lup. After imposing the 
absorbing column density found for IM Lup both spectra 
show the same features. The remaining small differences 
in the line fluxes are due to small abundance effects and 
IM Lup is a little stronger in the short-wavelength con- 
tinuum, because of its strong, hot corona. Notable differ- 
ences turn up, when the same experiment is done with the 
spectrum of TW Hya, because it lacks the strong emission 
component at 2.1 keV. 



5. Discussion 

5.1. X-ray properties 

IM Lup shows many characteristics of an active star. It has 
a bright, hot corona and its abundances follow the IFIP pat- 
tern. These traits are common to both CTTS and WTTS. 
In the density diagnostic of the He-like triplets there is no 
deviation from the coronal limit, however, due to the low 
signal a high-density state is allowed at the 90% confidence 
level. The signal of the O VII triplet is too weak to check for 
any soft X-ray excess. Thus, we conclude that - from the 
X-ray point-of-view - it is more likely that IM Lup shares 
its characteristics with WTTS, despite the presence of its 
disk. If confirmed, this suggests that the distinctive charac- 
teristics of CTTS are really due to accretion: An excess of 
soft X-ray emission originates in the high-density environ- 
mcnt of an accretion spot as simulated by iLamzinl (J1998) 
and IGunther et aLl (|2007h . Our observation is incompati- 
ble with any hypothetical scenario where the high-density 
signatures of the He-like triplets are attributed to the ac- 
cretion disk itself because they are not observed (on the la 
confidence level) in IM Lup, which has a disk. 

The only peculiarity in comparison to the WTTS sam- 
ples of COUP and XEST is the high bolometric luminosity 
of IM Lup. It belongs to the brightest objects of its kind. 
This may be due to an overestimation of the distance. We 
use 190 ± 27 pc from lWichmann et al.l (|l998f) based on the 
HIPPARCOS paralla x, but other metho ds lead to lower 
values of 140 ± 20 pc (jHughes et al.lll993l ). Both values are 
compatible within the 2a errors, but if the lower value was 
the true one, X-ray emission measure and bolometric lumi- 
nosity would be only half of the values given above. This 



does not change any of the other conclusions of our analysis, 
but IM Lup would fit in better with the bulk of the COUP 
and XEST objects in this case. Still, IM Lup is an active 
star with a bright and hot corona and log(L x / ' Lboi) = —3.1, 
similar to saturated MS stars and also typical for CTTS 
dGiidel fc Nazell200l . 

5.2. Mass accretion? 

IM Lup has a high level of chromospheric activity, which 
causes the emission observed in Ha. Other M dwarfs 
with active chromosphere also show H a emission, typically 
with an equivalent w idth of a few A (jWorden et al.lll98ll : 
iRobinson et al.|[l990l) . very similar to IM Lup. Their Ca H 
and K lines often show a very narrow emission peak, too. 
The Ha line profiles vary, but in most active stars there 
is a single, narrow emission peak, sometimes with a re- 
versal due to self-absorption in the centre. This is in line 
wi th the central peak i n the IM Lup data. In the s amples 
of IWordenet all (|1981[ ) and IRobinson et~aTl (|1990l ) there 
is not a single case of Ha in emission, where simultane- 
ously broad absorption is present. The same is true for 
WTTS , which have na rrow emission lines of symmetric 
shape (jFanget al.ll2009l ). 

In contrast, red-shifted absorption dips and red- and 
blue-shifted emission components a ppear in accretion mod - 
els in most viewing geometries (|Muzerolle et al.l |2001| ). 
Thus, in IM Lup the Ha line profile must be caused by addi- 
tional material which is not present in MS stars: Candidates 
are a circumstellar envelope, the disk or accretion funnels. 
We do not detect veiling of the photospheric emission lines, 
thus any mass accretion must be very low, if present at all. 
Still, the narrow emission components might originate in a 
weak accretion funnel. They are shifted by ±100 km s _1 , 
because we see accretion towards both, the observers side of 
the star and the back side. The line shift is smaller than usu- 
ally seen in emission lines from accretion funnels of CTTS, 
but given the uncertain stellar parameters of IM Lup this 
could be due to a combination of low mass, i.e. smaller 
gravitational potential, and viewing geometry, where the 
accretion funnel is not parallel to the line-of-sight. 

In some stars the Ha line profile is also influenced by 
strong flares but this cannot be the case in IM Lup, because 
no flare shows u p in the X-ray l ightcu rve. 

The study of lWhite fc Basril (|2003l) correlated the width 
of the Ha line at 10 % of the maximum with accretion. They 
suggest a full width at 10 % of the maximum > 270 km s _1 
as borderline between accreting and non-accreting sources, 
a slightly lower limi t of > 200 km s _1 was presented by 
Uavawardhana et al.l (|2003l) . The line width of IM Lup falls 
just between those two limits; it would be classified as non- 
accr eting by the one a nd as accreting by the other bound- 
ary. |Natti~eLal. (2004) showed that a quantitative relation 
between line width and mass accretion rate exists, which is 
valid from CTTS to sub-stellar objects. This relation yields 
an accretion rate of 10 _11 Mq yr _1 , but IM Lup is chro- 
mospherically very active, which explains a large fraction 
of the Ha flux. Additionally, not all relations employed to 
convert the Ha flux to a mass accretion rate agree, if ac- 
cretion is present at all, therefore we treat this number as 
a rough estimate only. 

In summary, the Ha EW classifies IM Lup as a WTTS, 
but the Ha line profile is complex as in CTTS. This sug- 
gests that mass accretion occurs at a very small rate. The 



H. M. Gunther et al.: The disk-bearing young star IM Lup 



energy released from accretion is too low to influence the 
X-ray properties significantly. For a mass accretion rate of 
10 —11 Mq yr -1 , at most 0.5% of the X-rays can be due 
to accretion, assum ing that 10% of the accretion energy is 
released in X-rays (jGiinther et al.l 120071 ) . For the interpre- 
tation of the X-ray data we can thus safely treat IM Lup 
as a non-accreting source. 



5.3. Varying absorption? 

The HARPS Ha line profiles are presented in Fig. [6] While 
the fitting with Gaussians (table [5]) yields an absorption 
component reaching over the entire line profile, for non- 
Gaussian components absorption is only required on the 
red side where the observed line profile falls below the con- 
tinuum. The situation is different in the ANU data, which 
was obtained simultaneous to the Chandra observation. No 
absorption is visible here on the first day and wide absorb- 
tion wings are present on both sides of the line on the second 
day. We postulate that matter from the disk on its way onto 
the star causes the Ha absorption. With IM Lup's mass and 
radius (Sect. [5]) we calculate the Keplerian rotation of the 
disk. At the co-rotation radius of 6 i?* it matches the stellar- 
rotation period from Sect. 14.41 in this region the inner disk 
should be truncated. The absorption component is strong 
at 100-400 km s _1 , which corresponds to the free-fall ve- 
locity from the co-rotation radius. Thus the Ha absorbing 
gas may be found in the gap between star and disk; how- 
ever, this scenario can only explain the red-shifted part of 
the absorption, but it strengthens the argument in Sect. 15.21 
that the Ha line width at 10% maximum could be inter- 
preted as accretion. We speculate, that the blue-shifted ab- 
sorption in the ANU data is caused by an outflow. CTTS 
commonly drive winds or even jets, where t he mass loss 
is of the order of 10% o f the accretion rate (jCabrit et al.l 
Il990t ICoffev et al.ll2008t) . The blue-shifted absorption ap- 
pears only on the second day of the ANU observation, when 
the absoprtion on the red side, which we interpret as accre- 
tion, is stronger. 

We divided the Chandra observation in two parts and 
fitted the rin indenpcdcntly to test if the absorption in the 
X-ray data changes with time. We found no change of the 
fitted value. Also, the hardness ratio of the X-ray spectra 
does not change during the observations. This does not rule 
out variability in n#, but is sets an upper limit. We esti- 
mate that an additional absorbing column of > 10 21 cm~ 2 
would have been detected. Does the accretion scenario ex- 
plain why strong Ha absorption coincides with weak X-ray 
emission? Given the weak limit on variability in nn and 
the fact that we have only two datapoints of simultaneous 
X-ray and optical spectra, the relation between Ha absorp- 
tion and X-ray luminosity is -at best- suggestive. Still it is 
intriguing to interpret the lower X-ray luminosity and the 
Ha line profiles as absoption by an accretion funnel. This 
should cause a hardening of the X-ray absorption. This is 
not observed, but the limit on extra absorption is weak due 
to the low count rates. 



allows the disk to reach closer to the star, i.e., weaker emit- 
ters accrete and become CTTS. 



5.4. IM Lup in context 



If taken seriously, the relation between X-ray luminosity 
and Ha absorption fits to the fact that CTT S are on average 
less luminous in X-rays than WTTS (e.g . iPreibisch et al.l 
20051: iTelleschi et al.ll2007bt ). iDrake etal . (2009) suggested 
that it is not the accreting nature of the CTTS that causes 
their low luminosity; rather, it is the low luminosity that 



ISchisano et al.l (|2009f ) observed variability in optical spec- 
tra of T Cha. Like IM Lup T Cha is also a transition ob- 
ject from CTTS to WTTS without veiling of photospheric 
lines, i.e. with little or no mass accretion. ISchisano et al.l 
(2009) find the Ha profiles to change from emission lines 
to absorption components over time, simultaneously they 
observe changes in the reddening and develop a scenario 
based on inhomogeneous circumstellar extinction similar to 
our ideas for IM Lup. They argue this to be caused by the 
inner parts of the disk, where grain growth takes place and 
possibly planets form. 

Only a small number of CTTS have been observed with 
X-ray grating spectroscopy, most of them show low f/i 
ratios in the He-like triplets, but T Tau does not. Does 
IM Lup confirm a CTTS population with low densities like 
T Tau or does it behave as WTTS? Although the complex 
Ha line profile points to mass accretion, IM Lup differs 
from T Tau in other characteristics. The main component 
of the multiple system T Tau is much more massive, close to 
the boundary to the HerbigAe/Be stars and it accretes at 
much higher rates. Whereas IM Lup's X-ray emission must 
be coronal, T Tau could power a significant part of its lu- 
minosity from accretion. We therefore expect IM Lup to 
resemble the WTTS and young MS stars in its X-ray prop- 
erties. Again, only very few WTTS are observed with X-ray 
gratings, but they all look similar to MS stars. Despite the 
small number, it is credible that their emission is driven 
exclusively by a solar-type corona and we find no evidence 
that the transition object IM Lup, which still posseses a 
disk, is different. The transition from a CTTS to a WTTS 
in X-rays seems to coincide with the decline of accretion as 
expected from accretion shock models for CTTS. However, 
the Ne IX triplet is the strongest discriminator and the orig- 
inal driver of the reported observation, but at the 90% con- 
fidence level the CTTS-like high-density scenario cannot be 
excluded for IM Lup. 



6. Summary 

We present a deep Chandra /HETGS grating spectrum of 
IM Lup, a WTTS with a circumstellar disk, and accom- 
panying optical spectroscopy, where some of the spectra 
were obtained simultaneous to the X-ray spectra. IM Lup 
is an active star with an IFIP effect in a hot corona. The 
Ne IX He-like triplet is in the low-density limit, but due to 
the low signal we also cannot exclude a low f/i scenario. 
It seems, that the X-ray signatures of CTTS, a soft ex- 
cess and high-density He-like triplets, are not correlated 
with the presence of the disk, but with the active accre- 
tion of CTTS. IM Lup's optical spectra exhibit emission 
lines, which are typical for chromospherically active stars. 
This fits well with the high X-ray luminosity of IM Lup. 
However, the Ha line profiles are more complex than ex- 
pected in MS chromospheres. This could possibly be due 
to weak accretion funnels as in CTTS with a low accretion 
rate. 
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